Study of the ground-state energy of 40 Ca with the CD-Bonn nucleon-nucleon potential 
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We have calculated the ground-state energy of the doubly-magic nucleus 40 Ca within the frame- 
work of the Goldstone expansion using the CD-Bonn nucleon-nucleon potential. The short-range 
repulsion of this potential has been renormalized by integrating out its high-momentum components 
so as to derive a low-momentum potential Vi ow -k defined up to a cutoff momentum A. A simple 
criterion has been employed to establish a connection between this cutoff momentum and the size of 
the two-nucleon model space in the harmonic oscillator basis. This model-space truncation approach 
provides a reliable way to renormalize the free nucleon-nucleon potential preserving its many-body 
physics. The role of the 3p — 3h and 4p — 4h excitations in the description of the ground state of 
40 Ca is discussed. 
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A fundamental goal of nuclear physics is to describe 
the properties of nuclei starting from the forces among 
nucleons. To this end, one has to employ many-body 
methods well suited to handle the strong short-range 
correlations that are induced in nuclei by the free-space 
nucleon-nucleon (NN) potential Vnn- In other words, 
the method employed should produce results which are 
only slightly affected by the approximations involved, in 
the sense that it should be possible to keep the latter 
under control by way of convergence checks. This first- 
principle approach to the description of nuclear structure 
is nowadays referred to as ab initio approach [l[ . 

In the last decade, thanks also to the considerable 
increase in computer power, a substantial progress has 
been achieved in microscopic approaches to the nuclear 
many-body problem, such as the Green's function Monte 
Carlo (GFMC) @, no-core shell model (NCSM) @, ||, 
and coupled-cluster methods (CCM) 

Historically, the first calculations on light p-shell nu- 
clei using the GFMC method were performed in the mid- 
1990s by the Argonne group Q. They employed the high- 
precision NN potential Q plus a three-body in- 
teraction, the latter being fitted to reproduce the binding 
energy of selected few-body systems. Since then, these 
calculations have been successfully performed up to mass 
A = 12 d, HI, which is the present limit of GFMC with 
the available computer technology, owing to the exponen- 
tial growth of the spin-isospin vector size (lfj| . This limit 
may be overcome by introducing effective interactions. In 
this way, the NCSM and CCM allow to perform calcula- 
tions beyond the p-shell mass region. More precisely, the 
NCSM has been applied 0,0] to nuclei with mass A < 16, 
using either coordinate-space or momentum dependent 
Vnn's 11]. In some cases, also three-body interactions 

have been included [H EE El- The CCM can Poten- 
tially be used for much heavier systems; in fact, in the 
late 1970s it was applied to the doubly-closed nuclei 16 
and 40 Ca [l5j]. Actually, coupled-cluster calculations em- 



ploying modern Vnn's |6j have been recently performed 
for valence systems around 16 O [H, El El- CCM-like 
calculations are also those performed for nuclei up to 16 O 
in Refs. [ljl[I(|, where different realistic NN potentials 
have been used in the framework of the unitary model- 
operator approach (UMOA) @, [U, H3] . 

In Refs. [23L [1H a method to renormalize the bare 
NN interaction has been introduced, which may be con- 
sidered an advantageous alternative to the use of the 
Brueckner G matrix. A low-momentum model space is 
defined up to a cutoff momentum A, and an effective 
low-momentum potential Vj ow _k, which satisfies the de- 
coupling condition between the low- and high-momentum 
spaces, is derived from the original Vnn- The V[ ow -k 
is a smooth potential which preserves exactly the low- 
momentum on-shell properties of the original Vnn and 
can be used directly in nuclear structure calculations. 

Recently, we have investigated [III [26| how A is re- 
lated to the dimension of the configuration space in the 
coordinate representation where our calculations are per- 
formed. We have introduced a simple criterion to map 
out the model space made up by the two-nucleon states 
in the harmonic-oscillator (HO) basis according to the 
value of the cutoff momentum A. The validity of this 
procedure was tested by calculating, in the framework of 
the Goldstone expansion, the ground-state (g.s.) energy 
of 4 He with the CD-Bonn [H, N 3 LO H, and Bonn A 
[291 ] potentials, and comparing the results with those ob- 
tained using the Faddeev-Yakubovsky method. Taking 
into account perturbative contributions up to fourth or- 
der in V[ ow _k, we have found that the energy differences 
are at most 390 keV. The limited size of the discrepancies 
shows that this approach provides a reliable way to renor- 
malize the NN potential preserving the physics beyond 
the two-body system too. 

We also performed calculations for heavier systems, 
such as 16 and 40 Ca, and obtained converged results 
for the CD-Bonn NN potential using a limited number 
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of oscillator quanta. As regards 40 Ca, the g.s. energy 
was calculated including Goldstone diagrams only up to 
third order in Mow-k- 
it seems fair to say that, at present, an ab initio cal- 
culation for 40 Ca represents a major step on the way to 
the fully microscopic description of nuclear systems be- 
yond 16 0. In this work, we improve our calculation of 
the 40 Ca g.s. energy including all the fourth-order con- 
tributions. A main motivation for this extension of our 
calculation of the g.s. energy of 40 Ca is to study the role 
of a higher-order class of excitations, namely the 3p — 3h 
and Ap — 4h ones, which come into play starting from the 
fourth order of the Goldstone expansion. It should be 
pointed out that this is the first fully microscopic study 
of this nucleus, apart from a very preliminary calcula- 
tion by Kumagai et al. [30( and Fujii et al. [20J in the 
framework of the unitary model-operator approach. 

The first step of our calculation is to renormalize the 
short-range repulsion of the NN potential by integrat- 
ing out its high momentum components through the so- 
called Viow-k approach (see Refs. 0, HH). The low- 
momentum potential Viow-k preserves the physics of the 
two-nucleon system, and consequently the 'x 2 j datum of 
the original NN potential, up to the cutoff momentum 
A. As mentioned above, it is a smooth potential and 
can therefore be used directly in a perturbative nuclear 
structure calculation. 

The Viow-k is defined in the momentum space, and it 
is desirable to relate the cutoff momentum A to the size 
of the harmonic-oscillator (HO) space in the coordinate 
representation [2^, [2(|, where we perform our calcula- 
tions for finite nuclei. 

Let us consider the relative motion of two nucleons 
in a HO well in the momentum representation. For a 
given maximum relative momentum A, the corresponding 
maximum value of the energy is: 



+ 2n b + l b < A max . (4) 

Making use of the above approach, in this paper we 
have studied the g.s. energy of 40 Ca within the frame- 
work of the Goldstone expansion [32j . We start from the 
purely intrinsic hamiltonian 

x ' i=l i<j 

where Vy stands for the renormalized Vat at potential plus 
the Coulomb force, and construct the Hartree-Fock (HF) 
basis expanding the HF single particle (SP) states in 
terms of HO wave functions. The next step is to sum up 
all the Goldstone linked diagrams for the ground-state 
energy up to fourth-order in the two-body interaction. 
The complete list of the fourth-order diagrams can be 
found in Ref. [H, [H| . Using Pade approximants [H, HH 
one may obtain a value to which the perturbation series 
should converge. In this work, we report results obtained 
using the Pade approximant [2|2], whose explicit expres- 
sion is 



Eq(1 + 7i + 72) + -Ei(l + 72) + E 2 



where 



E 2 Ej - E 2 E 3 + giTi 

^ 1 TP J? T?2 ' 771 



-C'max — -i 

where M is the nucleon mass. 

This relation may be rewritten in terms of the maxi- 
mum number -/V max of HO quanta in the relative coordi- 
nate system. For a given HO parameter huo we have: 

/ 3\ h 2 A 2 

(^V max + -J^ = — . (2) 

The above equation provides a simple criterion to map 
out the two-nucleon HO model space. If we write the 
two-nucleon states as the product of HO wave functions 

\a b) = \nJaja, n b l b j b ) , (3) 

our HO model space is defined as spanned by those two- 
nucleon states that satisfy the constraint 



Ei being the ith order energy contribution in the Gold- 
stone expansion. 

In principle, our results should not depend on the HO 
parameter fku, whose value characterizes the HO wave- 
functions employed to expand the HF SP states. Ac- 
tually, our calculations are made in a truncated model 
space, whose size is related to the values of the cutoff 
momentum A and the Tilo parameter by relations (|5J) and 
(ID). Obviously, for A max — > 00 this dependence disap- 
pears. So, we perform our calculations increasing the 
A max value (and consequently A) for different frui values. 
Finally, we choose the results which correspond to the 
value of huj for which they are practically independent of 
A max . 

In Fig. Q] we report the ground-state energy of 40 Ca 
calculated with the CD-Bonn potential [27j . The straight 
red line indicates the experimental datum [36| while the 
other curves represent our calculated values, for differ- 
ent values of hu>, versus the maximum number of HO 
quanta A max that limits the two-nucleon configurations 
according to relation Q. 
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FIG. 1: (Color online) Ground-state energy of 40 Ca with the 
CD-Bonn potential as a function of iV max , for different values 
of hui. The straight line represents the experimental value, 
while the dashed one our converged result. The value of the 
energy difference between the former and the latter is also 
reported. 



We obtain convergence for Hlo = 25.5 MeV, and the 
corresponding energy is (— 314 ± 3) MeV, as indicated in 
Fig. [T]by the dashed line. The error has been evaluated 
taking into account the dependence of our results on h\o 
[26| . In Ref. [25[ we calculated the 40 Ca g.s. energy tak- 
ing into account all the contributions in the Goldstone 
expansion up to third order. The converged value, ob- 
tained with the Pade approximant [2|1] turned out to be 
(-308 ± 3) MeV with huj = 25.5 MeV. The difference 
between the two results is only about 2% of the total 
binding energy, the converged value at fourth order be- 
ing 6 MeV more attractive than the third-order one. This 
is in line with the outcome of our calculations (2f| H(| of 
the g.s. energy of 4 He and 16 with the CD-Bonn poten- 
tial. In those cases the converged fourth-order result was 
0.6 MeV and 3 MeV more attractive, respectively, than 
the converged third-order one. 



TABLE I: Calculated 2p-2h, 3p-3h, and Ap-Ah fourth-order 
contributions (in MeV) to the g.s. energy of 40 Ca with the 
CD-Bonn potential. Calculations correspond to jV ma x = 10. 



Nucleus 


2p - 2h 


3p — 3h 


Ap — Ah 


4th order 


40 Ca 


-15 


-16 


+24 


-7 



It is worth now to make a brief discussion about the 
role of 3p— 3h and Ap—Ah contributions which, as pointed 
out before, come into play only at the fourth order of the 
Goldstone expansion and beyond. To this end, we report 
in Table [J the 2p — 2h, 3p — 3ft., and Ap — Ah fourth-order 
contributions to the g.s. energy of 40 Ca with the CD- 
Bonn potential. Since we use the Hartree-Fock basis, all 
Xp — lh excitations of the ground state are identically 
zero. The inspection of Table Q] evidences that the role 
played by the 3p— 3h and Ap—Ah excitations is significant, 
their net repulsive contribution being 8 MeV. The latter 
counterbalances the fourth-order contribution of the 2p — 
2h excitations, so that the total fourth-order contribution 
is only -7 MeV. 

As regards the comparison with experiment, our cal- 
culated 40 Ca binding energy underestimates the experi- 
mental one by 28 MeV, which is 8% of the experimental 
binding energy. Our calculations of 4 He and 16 O with the 
CD-Bonn potential [25|, [2|| show the same percent differ- 
ence from the experimental data. This seems to confirm 
the need of a three-body force in addition to the NN 
CD-Bonn potential, in order to compensate for the lack 
of attraction of the latter. 

In conclusion, this work presents a fully microscopic 
calculation of the ground-state energy of 40 Ca with the 
CD-Bonn NN potential. We hope that this may stim- 
ulate, and provide some useful hints to, future non- 
perturbative ab initio calculations. 
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